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Abstract 

We have analyzed the morphology of ferroelectric domains in very thin films of multiferroic 
BiFeOa. Unlike the more common stripe domains observed in thicker films BiFeOs or in other 
ferroics, the domains tend not to be straight, but irregular in shape, with significant domain wall 
roughening leading to a fractal dimensionality. Also contrary to what is usually observed in other 
ferroics, the domain size appears not to scale as the square root of the film thickness. A model is 
proposed in which the observed domain size as a function of film thickness can be directly linked 
to the fractal dimension of the domains. 
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Multiferroic materials are currently attracting a great deal of attention on account of 
their interesting physics and potential applications Of these, possibly the most studied 
is the perovskite BiFeC>3 (BFO) , one of the very few materials which is multiferroic (ferro- 
electric and antiferromagnetic) at room temperature j^, Q|. Its lead- free nature and large 
remanent polarization have already motivated Fujitsu to use it a S the active layer in pro- 
totype ferroelectric memories [5fl; also, sub-lattice magnetic switching using voltage has been 
demonstrated jf]], which may find its way into spintronic applications via exchange bias 
The possible coupling between ferroelectric and antiferromagnetic domains has triggered a 
flurry of work on the morphology and functional properties of the domains jf], [s[ fol.Tlol. 
The ferroelectric domains are generally found to be strai ght- walled, and to follow the well- 
known scaling law of Landau, Lifshitz and Kittel (LLK) [I2I I13I. fl4 ]. that is, their domain 
size grows proportionally to the square root of film thickness [1 Qt ] . 

BFO has a complex phase sequence, both in bulk and in thin film form. In bulk, it goes (on 
cooling) from cubic to orthorhombic to rhombohedral, with a recently-discovered possible 



metal-insulator transition at high-temperature [15|. The room-temperature rhombohedral 
phase is normally monoclinic for epitaxial thin films, but it appears to be tetragonal below 
a critical thickness, which for BFO grown epitaxially on SrTiOs substrates is of the order 
of lOOnm . We have analyzed the morphology and scaling of the domains in the small- 
thickness (tetragonal) regime and found it to be fundamentally different from that observed 
at higher thickness: i) the domains are not straight, but irregular in shape; ii) there is 
significant domain wall roughness with a fractal dimensionality; and iii) the average domain 
size does not scale as the square root of the film thickness. Taking into account the fractal 
nature of the domain walls, we have derived a new relation which accounts for the anomalous 
domain size behaviour, with the scaling exponent being a function of the fractal Hausdorff 
dimension. 

The thin films of BFO were grown by pulsed laser deposition on SrTiOs with a conductive 
buffer layer of (La,Sr)Mn03 [Hfl, and the ferroelectric domain morphology was studied by 
way of piezo-response atomic force microscopy (PFM). In Figure 1 we show the surface 
topology and PFM response of some films. The shape of these "virgin" domains is highly 
irregular, and not correlated with topo grap hy. Irregularly-shaped domains have also been 



seen in much thinner films (30nm or less 



Ufl), so it appears to be a general property of ultra- 



thin BFO, with the critical thickness depending on the choice of substrate and electrode, as 
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FIG. 1: Surface topography and PFM scans of domains in BFO films of different thickness. The 
PFM cantilever was oriented along a (100) axis 

well as growth conditions. 

Figure 2 shows the average domain period as a function of film thickness. For comparison, 
equivalent data for domains in other ferroic systems (ferroelectric or ferromagnetic) are 
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While ferroelectric domains are generally smaller than ferromagnetic domains 



2l| . the ferroelectric domains in multiferroic BFO are noticeably bigger than 
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those in "pure" ferroelectrics. This suggests a higher energy cost of the domain walls 
22]], possibly due to the contribution from the magnetoelectric coupling term. Even more 
surprising is the fact that the domain periodicity does not scale as the square root of the 
film thickness. In our case, the exponent for the thickness dependence in w = Ad 1 (where w 
is domain period and d is film thickness) is around 7 ~ 3/4. This is not due to the fact that 
the domains are of mosaic-type rather than stripes, since such domains are also predicted 



to follow the 7=1/2 scaling [l4j, so a different explanation is required. 



In order to gain further insight into this problem we have studied the morphology of 
the domain walls. Firstly, using the program WSxM H|, we have analyzed the ratio of 



perimeters (P) to areas (A) for the spontaneous domains [24]. The scaling is such that 
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FIG. 2: Above, periodicity of in-plane and out-of-plane domain contrast in the PFM scans as a 
function of BFO film thickness. The straight lines are least-fit squares giving scaling exponents 
of 0.70±0.07 (out-of-plane) and 0.81±0.04 (in-plane). Below, comparison between the ferroelectric 
domain periodicity of BFO and the domain periodicity of other ferroics found in the literature. 

P oc A H ^ 2 , where H\\ is the Hausdorff dimension of the domain walls in their projection 
on the x-y plane accessible by the PFM; if the domain walls were perfectly smooth, ii"||=l, 
otherwise 1 < H\\ < 2. Plotting log(P) vs \og(A) (figure 3), H\\ = 1.45 ± 0.2 was obtained. 
This value is in close agreement with that found in BFO-doped lead zirconium-titanate 
(PZT) films, where fractal ferroelectric domains were written using a PFM tip via liquid 



electrodes 
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To relate these observations to the nature of the domain walls, we have also directly 
measured the domain wall roughness of AFM-written linear domains, created by alternate 



applications of negative and positive voltage (± 8V) to the film surface. The resulting ( 
walls were analyzed using the pair-correlation method described by Paruch et al 



omam 



26|]. The 
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correlation function of relative displacements essentially measures the local variance of the 
domain wall position from an elastically ideal flat configuration as a function of the length 
L along the wall, and is predicted to show a power-law growth at equilibrium, governed by a 
characteristic roughness exponent zeta. We observe such a power law growth at short length 
scales (L<100 nm) followed by a saturation of the correlation function. The value of the 



roughness exponent ( ~ 0.5 — 0.6 is higher than that observed in PZT [26j, and could be an 
indication of either a lower dimensionality (being close to the theoretically predicted value 
of C=2/3 for a one-dimensional elastic domain wall in random bond disorder potential), or 
possibly a different type of disorder playing an important role in a system with combined 



magnetic and electric effects 



27|]. 



There are additional self-consistency checks possible for the fractal Hausdorff dimension 



requiring time-resolved data for current transients 



281 ] . For well-annealed PZT, analysis of 



I(t) yielded [29|] H = 2.55 ± 0.25, compatible with the present results on BFO assuming 
straight walls normal to the film surface, that is, H = H\\ + H±=l. 5+1.0 where H\\ and 
H± are, respectively, the in-plane and out-of-plane components of the fractal dimension. A 
similar result was also inferred for SBN (strontium barium niobate) by Kleemann et al. J30I ]. 
who observed a dimensionality of 1 < H—l < 1.7, inferring H = 2.35±0.35. The dimension- 
ality of 2.5 is also compatible with hyperscaling near critical points in ferroelectrics j^ll, 



Finally, it is also worth mentioning the work of Yu and Randall [33|, who showed fractal-like 
(dentritic) domain structures in relaxor ferroelectrics, although the fractal dimension was 
not calculated in that case. 

The fractal nature of the domain walls, and even the specific Hausdorff dimension of ~ 2.5, 
appear therefore to be rather more general than previously thought, so it is important to 
understand how fractal walls may affect domain scaling. A fractal distribution of domain 
sizes (Cantor sets) has been analysed [34] for ferroelectric stripe domains, but the walls in 
these studies were straight, and a conventional LLK exponent 7 = 1/2 was still found, unlike 
that of our films. Here, instead, we have modified the domain size scaling problem in order 
to incorporate the fractal roughening. The LLK law arises from the need to minimize the 
total energy of the domains and that of the domain walls. The energy density of the domains 
is proportional to their size w (irrespective of whether they are stripe-type or mosaic-type 

m 
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E domain = Uw (l) 

Where U is a constant arising from either depolarization, demagnetization, strain, or 
a combination of them. The energy of the domain walls is proportional to their number 
density, which is found dividing the area of the film by the area of each domain, i.e., N = 
L 2 /u> 2 =number of domain walls. The energy of each domain wall is equal to its energy 
density (a) times the area of the wall, which is essentially the domain perimeter times the 
thickness of the film: aPd. If the domain walls are fractal, then the domain perimeter does 
not scale linearly with the average domain size w, but fractally as P = w Hn , where H\\ is 
the Hausdorff dimension. Hence E wa u = aw H ^d. The energy of all the domain walls in the 
crystal is thus E wa u x N = crw H ^dL 2 /w 2 , and the energy density (per unit area of film) is 
found dividing by L 2 . Thus: 

E wa lls = 0"W H U A; (2) 

ur 

The optimum domain size is found by adding (1) and (2) and minimizing with respect to 
w. This leads to 



w 



(2 - H\\) ^ ' d\ z ' H w ) = kd\ 3 ~ H W J (3) 



The above treatment assumes that the domain walls are straight in the vertical direction. 
Otherwise, the vertical size of the domain walls would scale as d H± , and Eq. [3] would become 
instead: 



w = kd\ 3 - H wJ (4) 

For H±=H»=1 (flat walls) we recover the standard Kittel result 7 = 1/2. 

PFM does not allow establishing, whether the walls are rough or straight along the z- 
direction. Previous work in PZT [26] and Bi 4 Ti 3 0i2 j35] shows that domain walls may be 
rough along one direction and flat along another, due to anisotropy in the strength of the 
dipole-dipole and elastic interactions. Here, comparison between eq. El eq. 31 and the 
measured domain periodicity as a function of thickness can help to indirectly determine H±. 

For the measured H\\ ~ 3/2, assuming straight walls along the vertical direction, eq. [3] 
predicts a value 7 ~2/3 for the exponent in w = kd 1 . If, on the other hand we assume that 
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FIG. 3: Above, perimeter as a function of area for spontaneously formed domains. The slope of 
the log-log plot is related to the Hausdorff dimension: P oc A H ^ 2 . Below, correlation function for 
the roughness of the written domain walls as a function of lengthscale, measured for a 70 nm thick 
BFO film. An average value of the roughness exponent exponent £=0.56 is obtained for the seven 
domain walls indicated. This higher value of £ for BFO compared to that for PZT suggests either 
a lower dimensionality in a random bond pinning scenario, or the presence of stronger individual 
pinning centers which could also be responsible for the fractal structure of the spontaneous domains 
in the BFO films. 

H«=H±, then eq. [4] predicts 7=1 instead (although it is worth mentioning that values of 
7 approaching 1 would also be consistent with the classical theory for ferroelastic domains, 



which predicts 7=1 for small thickness 



36|). The actual measured value for the periodicity 



of the in-plane domains is (see Figure 2) 7 RiO.7-0.8, which suggests that the domain walls 
are anisotropically roughened, being more irregular in the horizontal direction than in the 
vertical one. 

The present theoretical framework is valid not just for BFO, but for any ferroic (ferro- 
electric, ferromagnetic or ferroelastic) material with fractal walls. Although it is consistent 
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with the empirical results, the domains in the thinnest films are close to the resolution limit 
of PFM, and thus there is a substantial error bar associated with their measured periodicity. 
Also, the value of the scaling exponent 7 appears to be different for out-of-plane and in-plane 
variants (Figure 2-a), although the difference is within the statistical error. Finally, while 
the fractal analysis was performed on the samples with best signal-to-noise ratio, there may 
still be some electronic noise contribution to the fractal dimension. 

In summary, we have analyzed the domain periodicity and domain morphology of very 
thin films of BFO and found: i) the spontaneous ferroelectric domains in BFO are generally 
bigger than those in other ferroelectrics, consistent with an increased energy cost of the 
domain walls due to magnetoelectric coupling; ii) the domains are irregularly shaped and 
characterized by a fractal dimension; iii) the domain size does not scale as the square root 
of the film thickness; and iv) the unusual scaling can be directly related to the the fractal 
dimension of the walls. Specific further tests on the fractal dimensionality inferred are 
suggested, particularly the use of switching current transients. We hope that these findings 
will motivate more research into the physics of domain walls in ferroic and multiferroic 
systems. 
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